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We consider the stochastic propagation of high-energy protons and nuclei in the cosmological microwave and 
infrared backgrounds, using revised photonuclear cross-sections and following primary and secondary nuclei 
in the full 2D nuclear chart. We confirm earlier results showing that the high-energy data can be fit with a pure 
proton extragalactic cosmic ray (EGCR) component if the source spectrum is oc E~^-^. In this case the ankle 
in the cosmic ray (CR) spectrum may be interpreted as a pair-production dip associated with the propagation. 
We show that when heavier nuclei are included in the source with a composition similar to that of Galactic 
cosmic -rays (GCRs), the pair-production dip is not present unless the proton fraction is higher than 85%. In 
the mixed composition case, the ankle recovers the past interpretation as the transition from GCRs to EGCRs 
and the highest energy data can be explained by a harder source spectrum oc E^^-^- E^^-^, reminiscent of 
relativistic shock acceleration predictions, and in good agreement with the GCR data at low-energy and holistic 
scenarios. While the expected cosmogenic neutrino fluxes at high energy are very similar for pure proton and 
mixed composition hypothesis, the two scenarii predict very different elongation rates from 10^^^ to 10^" eV. 



1. Introduction 

One of the keys to understanding the origin of cosmic -rays is the spectral shape of the transition from Galactic 
cosmic-rays to extragalactic cosmic rays. The Galactic origin of low-energy CRs is generally accepted, while 
the highest energy CRs (above ^ 10^^ eV) are no longer confined by Galactic magnetic fields and most 
probably originate from other galaxies. Thus, a transition between the two components has to occur in some 
energy range. The most natural location for this transition is around ^ 3 10^^ eV at a feature in the CR spectrum 
known as the ankle. This is the only energy range where the spectrum gets harder (i.e., its logarithmic slope 
gets smaller), offering a simple transition scheme to a harder CR component which is subdominant at lower 
energies. 

A different conclusion has recently been proposed on the basis of composition results from the HiRes Collab- 
oration, tentatively showing a transition from heavy to light primary nuclei at an energy around 5 lO^'^ eV 1 1 1, 
which can be identified with a second knee feature in the spectrum. Such composition measurements rely on 
statistical determinations of the elongation rate of extensive air-showers and depend on shower development 
simulations which are still quite uncertain. In this case the ankle is interpreted as an e+-e^ pair-production dip 
||2j resulting from proton propagation over large distances in the cosmic microwave background (CMB) and 
the transition from GCR to ECGR does not have to be at the ankle. For these second-knee-transition models 
the EGCR component is best fitted by assuming a power-law source spectrum in E^^^ with j3 ~ 2.6 |3 1 or 
/3 ~ 2.7 |2|. These studies, however, assume that EGCRs are only protons. 

In this study, we investigate how these conclusions are modified when nuclei heavier than protons are also 
present in the source. We apply the revised scheme of photo-nuclear interactions described in 0. Interest- 
ingly, we find that the two main aspects of EGCR phenomenology are jointly affected, namely the best fit 
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slope of the source spectrum and the interpretation of the ankle. The interpretation of the ankle as a pair pro- 
duction dip is not viable in the mixed composition case, unless the fraction of protons is above 85% . When 
reasonable assumptions for the fraction of nuclei in the EGCR are introduced, the traditional interpretation of 
a GCR/EGCR transition at the ankle is recovered and the best fit spectral index for EGCR is ^ 2.2 - 2.3, as 
expected from relativistic shock acceleration. Finally, the best fit EGCR spectrum occurs when both the GCR 
and EGCR component have the same spectral index, allowing for holistic approaches to the CR spectrum. 

2. Phenomenological source model 

Although extragalactic magnetic fields could modify the UHECR spectrum, we assume here that their intensity 
is low enough to be negligible and focus on a comparison between pure proton and mixed-composition models. 
The choice of a source composition is arbitrary in the absence of a source model, but we consider a "generic 
composition" assuming that EGCRs have the same relative source abundances as the best known, low-energy 
GCRs, derived from Ulysses and HEAO 3 data Q. We convert the differential abundance ratios, Xi, at a given 
energy per nucleon (E/A) into ratios at a given energy, = XiA^~^, where (3 is the spectral index, so that the 
GCR source spectrum is: Ni{E) oc ^iE^^. 

Since our goal is also to investigate the viability of the proton pair production dip as an explaination for the 
ankle when heavier nuclei are emited at the sources we will as well consider proton dominated compositions 
with small fraction of heavier nuclei. We will study in particular the impact of small fractions of Helium and 
Iron nuclei. 

Finally, we choose a rigidity dependent maximum energy at the source for the various nuclear species i?niax (^X) . 
We assume that energy losses and photo-fragmentation inside the source can be neglected, so that all nuclei 
with the same gyroradius behave the same way. This implies E'maxl^X) = Z x Ej^axiiH)- 

3. Results and discussion 

We compute the propagated spectrum of EGCRs for a uniform source distribution with negligible magnetic 
fields, using a Monte-Carlo technique. For this purpose we take into account all the relevant energy loss 
processes for protons and nuclei ^ In Fig.^, we show the resulting propagated spectra obtained with the 
prescription described above where E'maxdH) = 3 10^° eV. The error bars represent the l-cr fluctuations of 
the flux in each energy bin, as expected for a data set with the statistics of the AGASA experiment, i.e., 866 
events above 10^^ eV 0. We compare our results with the data from HiRes and AGASA experiments. In the 
pure proton case, we find again that the best fit spectrum has /3 ~ 2.6 and a pair production dip, as in previous 
work. Figures[3) and[3; show how the spectrum around the ankle and the GZK bump evolves with the fraction 
of protons in the source composition (above 10^'* eV) for two composition hypothesis. When the remaining 
composition is dominated by Helium, for fractions lower than 75%, the pair production dip is strongly atten- 
uated and the ankle is actually reversed for fractions lower than 60%. When small fraction of iron nuclei 
are present, the proton pair production dip is also modified but in a different way (due to the difference in 
the energy losses between He and Fe 1 1 11) the extragalactic spectrum appears flatter below 5 10^* eV and a 
galactic component is then needed to reproduce the low energy part of the ankle. Therefore, the interpretation 
of the ankle as an e+-e^ pair production dip requires a very large fraction of protons at the source. In contrast, 
Fig-Q shows two good fits to the highest energy data using mixed composition models with 13 = 2.2 or 2.3, 
corresponding to 50% or 40% protons at the source, respectively. 

^For further details about the treatment of these interactions and further references see f^. 
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Figure 1. Top left: Propagated spectra for sources with protons only and (3 — 2.6, and with mixed compositions and 
/3 = 2.2 and 2.3. Data points from the main experiments are indicated, with an arbitrary downward shift of AGASA's 
energy scale by 20%. Top right: evolution of the spectrum with the fraction of protons at the source, for an E~^'^' 
source spectrum and a heavier compenent dominated by Helium nuclei. Bottom left: Same as top right but with a heavier 
component dominated by iron nuclei. Bottom right: Elongation Rates corresponding to the pure extragalactic proton and 
mixed composition cases for the QGSJet and SIBYLL hadronic models, compared with Hires, Fly's eye and Yakutsk's 
results. 



The pure proton model with /3 = 2.6 provides a good fit to the data down to ^ 10^^ eV (assuming a uniform 
source distribution and no magnetic field), the transition from GCR to EGCR should occur before the ankle. 
However, shock acceleration processes are known to yield a spectrum with /3 ~ 2.2-2.3. As shown in Fig.^, 
such spectral indexes can provide equally good fits to the high energy data, assuming a realistic source model 
with a mixed composition. In such a scenario, the transition from GCRs to EGCRs occurs at the ankle, which 
thus keeps its "standard" interpretation. In addition, injection spectrum ~ E^^-^ has also been considered as 
the best fit to the low-energy data for the GCR component (e.g. 1 8 , 9 1). It is thus particularly interesting to note 
that assuming a similar composition for both the GCR and EGCR components is precisely what also makes a 
similar source spectrum possible. In addition, it has been shown in |7 1 that a source power-law index of ^ 2.3 
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is a necessary condition for holistic models, in which the same sources produce the CRs at all energies. The 
results presented here show that a self-consistent model can be built for CRs with a similar composition and 
spectrum at all energies provided the GCR/EGCR transition is at the ankle. 

Furthermore, we also calculated the expected cosmogenic neutrino fluxes for the pure proton and mixed model 
II II . For both cases we found very similar high energy neutrino peaks. This can be understood easily since 
protons dominate the high energy spectrum in the mixed case. We can then conclude for both composition 
hypothesis high energy neutrino fluxes should be detectable with experiments like Auger or Anita, but these 
fluxes could certainly not be used to distinguish between the two models. 

In Fig. [Ql, we show the predicted elongation rates for the composition hypothesis of Fig. [TJ, according to 
10000 Aires 1 10 1 simulated showers of the most abundant species II II (i.e, H, He, CNO, Mg, Si, Fe) between 
10^^-^ and 10^" eV. For this purpose we assume that the galactic iron component above 10^^-^ eV is the 
difference between the predicted extragalactic component and the experimental data. In the case of the pure 
proton extragalactic component, a very steep elongation rate is expected between 10^^-^ and 10^® eV in very 
good agreement with Hires data 1 1 1. As it is shown on Fig.^ and[n;, since the proton fraction as to be very 
high in the case of a second knee transition model (i.e a dip due to pair production), the elongation rate is 
expected to a have the characteristic shape of a rapid transition from a heavy galactic composition (shown by 
Kascade results | TS'l) to a very light extragalactic component unavoidably close to a pure proton expectation. 
Above 10^^ eV the expected elongation depend strongly on the hadronic models, the consistency of Hires 
data with this scenario is thus diffucult to establish, one can at least notice that the pure proton case is within 
the ~ 25 g.cm^'^ HiRes systematic error and the transition is in the expected energy range. For the mixed 
composition a continuous transition, from galactic iron to mixed composition (with a flattening between 10^®-^ 
and ~ 10^^-^ eV) and finally above 10^^-^ eV almost pure proton, is expected. This predicted shape is in good 
agreement (especially for the SIBYLL model) with Fly's eye and also Yakutsk results 11121 except below 10^* 
eV where the exact composition of the galactic component becomes relevant. 

Although the predictions are hadronic model dependant, the expected elongation rates appear extremely dif- 
ferent whatever the detail of the astrophysical hypothesis 111 II for the two transition models (i.e second knee 
or ankle). The current disagreements between different experimental results do not allow a conclusion at this 
point since each model can be favored or as well almost ruled out, whether one trusts Hires or Fly's eye elon- 
gation rate. However,the Pierre Auger Observatory should distinguish between the two hypothesis with high 
statistics and high accuracy hybrid elongation rate and finally explain the origin of the ankle. 
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